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INTRODUCTION 



A. Imaging the Nucleus 



We have a quantitative understanding of the strong interaction processes at extreme 
short distances in terms of perturbative QCD. We also understand the long distance prop- 
erties of hadronic interactions in terms of chiral perturbation theory. At the intermediate 
scale: the scale of quark confinement and the creation of [ordinary] mass we have an under- 
standing of numerous observables at about the 20% level from lattice QCD calculations and 
semi-phenomenological models. This is an extremely impressive intellectual achievement. 
However, the questions we have today about nuclear physics, are questions at the 1%, or 
even 0.1%, level relative to the confinement scale Aqcd ~ 300 MeV/c. For example, the n-p 
mass splitting of 1.3 MeV and the Deuteron binding energy of 2.2 MeV are < 1% effects that 
are crucial to the evolution of the universe. The patterns of binding energies of neutron and 
proton rich nuclei are even smaller effects, and are crucial to the synthesis of elements Z >Fe 
in supernovae and other extreme events. It is the QCD dynamics at the distance scale of 
1/Aqcd that gives rise to the origin of mass. To improve our understanding of confinement 
and of the origin of mass, we cannot rely solely on improvements in theory. We must have 
experimental observables of the fundamental degrees of freedom of QCD-the quarks and 
gluons-at the distance scale of confinement. The generalized parton distributions (GPDs) 
are precisely the necessary observables. 

Measurements of electro-weak form factors determine the spatial structure of charges 
and currents inside the nucleon. However, the resolution scale is not independent of 
the distance scale 1/a/Q^ probed. Deep inelastic scattering of leptons (DIS) and related 
inclusive high hadron scattering measure the distributions of quarks and gluons as a 
function of light cone momentum fractions, but integrated over spatial coordinates. Ji Q|, 
Radyushkin pl, and Miiller et al. 0|, defined a set of light cone matrix elements, now 
known as GPDs, which relate the spatial and momentum distributions of the partons. This 
allows the exciting possibility of determining spatial distributions of quarks and gluons in 
the nucleon function of their wavelength. 

The factorization proofs established that the GPDs are experimentally accessible 

through deeply virtual Compton scattering (DVCS) and its interference with the Bethe- 
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vcs 




Bethe-Heitler 



Figure 1: Lowest order QED diagrams for the process ep ep'y, including the DVCS and Bethe- 
Heitler (BH) amplitudes. The external momentum four-vectors are defined on the diagram. The 
virtual photon momenta are = [k — k'Y in the DVCS- and A'^ = {q — q'Y in the BH-amplitudes. 



Heitler (BH) process (Fig. [T]). In addition, the spatial resolution (Q^) of the reaction is 



independent of the distance scale ~ l/-\/tmin — t probed. Quark-Gluon operators are clas- 
sified by their twist: the dimension minus spin of each operator. The handbag amplitude 
of Fig. [21 is the lowest twist (twist-2) 7*p 'jp operator. The factorization proofs confirm 
the connection between DIS and DVCS. The proofs therefore suggest (but do not establish) 
that, just as in DIS, higher twist terms in DVCS will be only a small contribution to the 
cross sections at the and xbj range accessible with electrons from 6-12 GeV. 

In the formalism we are using [h|, the matrix elements of operators of twist-ra are 
independent (except for \n{Q^ / A'^^y)) evolution), but all observables of twist-n operators 
carry kinematic pre-factors that scale as [(t^in - t)/Q^f/^, [(-^)/(5^]"/^ or [M^/Q^f/'^. 
Diehl et al, [3| showed that the twist-2 and twist-3 DVCS-BH interference terms could be 
independently extracted from the azimuthal-dependence (0^7, §11 Ajl of the helicity depen- 
dent cross sections. Burkardt showed that the t-dependence of the GPDs at = is 
Fourier conjugate to the transverse spatial distribution of quarks in the infinite momentum 



extended 



frame as a function of momentum fraction. Ralston and Pire ISlI and Diehl 
this interpretation to the general case of ^ 7^ 0. Belitsky et al., [ll|| describe the general 
GPDs in terms of quark and gluon Wigner functions. 

These elegant theoretical concepts have stimulated an intense experimental effort in 
DVCS. The HI Q and ZEUS |l4| collaborations at HERA measured cross sections 
for xbj ~ 2,^ ^ 10~^. These data are integrated over and are therefore not sensitive to 



the BH-DVCS interference terms. The CLAS 
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and HERMES 
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collaborations mea- 
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Figure 2: Leading twist j*p — > amplitude in the DVCS limit. The initial and final quarks carry 
light-cone momentum fractions a; + ^ and a; — ^ of the light-cone momenta (1 ± -|- (Eq. 

The crossed diagram is also included in the full DVCS amplitude. The invariant momentum 
transfer squared to the proton is t = A^. 

sured relative beam helicity asymmetries. HERMES has also measured relative beam charge 
asymmetries and CLAS has measured longitudinal target relative asymmetries 

The HERA and HERMES results integrate over final state inelasticities of Mx < 2.9 GeV^ 
(or greater). Relative asymmetries are a ratio of cross section differences divided by a 
cross section sum. In general, these relative asymmetries contain BH-DVCS interference 
and DVCS^DVCS terms in both the numerator and denominator (the beam charge asym- 
metry removes the DVCS^DVCS terms only from the numerator). Absolute cross section 
measurements are necessary to obtain all DVCS observables. 

In Hall A, EOO-lloj2^ and E03-106[2l|, we measured absolute cross sections for H(e, e'7)p 
and D{e,e''j)pn at xbj = 0.36. The following section describes the methods and results of 
EOO-110. We propose to continue this program using the same experimental techniques, 
while taking advantage of the higher beam energy available with CEBAF at 12 GeV. 
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B. Review of Hall A DVCS EOO-110 Methods and Results 

The first draft pubhcation of EOO-110 can be found in Ref. j^l- Using a well understood 
experimental apparatus, this experiment measured both the unpolarized and polarized cross- 
section of the ep — epy process in the Deeply Virtual Compton Scattering (DVCS) regime 
at XBj = 0.36 and for at 1.5, 1.9, and 2.3 GeV^ . With electrons detected in the HRS- 
L, we have an absolute acceptance for electron detection understood to 3%, and a precise 
measurement of the scattering vertex and the direction of the virtual photon. DVCS is a 
three body final state, but at high and low t, the final photon is highly aligned with the 
virtual photon and therefore highly correlated with the scattered electron. Thus, even with 
a modest calorimeter, our coincidence acceptance for DVCS is essentially limited only by 
the electron spectrometer. As a consequence, very high values of the product of luminosity 
times coincidence acceptance are possible. The radiation hard PbF2 calorimeter gives a fast 
(Cerenkov) time response, and each channel is recorded with a 1 GHz digitizer which allows 
off-line identification of the DVCS photon. The identification of the exclusive channel is 
illustrated in Fig. [31 In Hall A, our systematic errors are minimized by the combination of 
the Compton polarimeter, the well-understood optics and acceptance of the High Resolution 
Spectrometer (HRS), and a compact, hermetic, calorimeter. All of those factors allowed the 
measurements of cross-sections. The high-precision electron detection minimizes systematic 
errors on t and (f)-y^. Therefore, we exploit the precision 0^^-dependence to extract the 
cross section terms which have the form of a finite Fourier series modulated by the electron 
propagators of the BH amplitude. 

The strong point of experiment EOO-lOO is that we measured absolute cross-sections. For 
example, the sin(0) term (modulated by BH propagators) of the helicity dependent cross 
section measures the interference of the imaginary part of the twist-2 DVCS amplitude with 
the BH amplitude, with a small contribution from an additional twist-3 bilinear DVCS term. 
The Q^-dependence of this term places a tight limit on the contribution of the higher twist 
terms to our extraction of the "handbag" amplitude. The unpolarized cross section is a sum 
of the BH cross section, the real part of the BH-DVCS interference, and a twist-2 bilinear 
DVCS term. The EOO-110 results show that the unpolarized cross section is not entirely 
dominated by the BH cross section, but also has a large contribution from DVCS. Explicit 
twist-3 terms were also extracted from the sin(20^^) and cos(20^^) dependence of the cross 
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Figure 3: Missing mass squared distribution for the H(e, e'7)X reaction, obtained in Hall A in 
EOO-110. The left plot shows the raw spectrum (after accidental subtraction) and the statistical 
estimate of the contribution of ep — > cYtt^ events of the type H(e, €'7)7^. This estimate was 



obtained from the ensemble of detected H(e, e'Tr^)Y events, with both photons from the vr' 



77 



decay detected in the calorimeter. The right plot shows the H(e, e'j)X spectrum, after the vr^ 
subtraction. 

sections. These contribute very little to the cross sections. 

In the kinematic regime of EOO-lOO, neither the extracted Twist-2 or Twist-3 observables 
show any statistically significant dependence on . This provides strong support to the 
original theoretical predictions that DVCS scaling is based on the same foundation as DIS 
scaling. We note that in the range 0.2 < xbj < 0.6, both the hi gher twist terms and 



24,^. Thus we 



InQ'^ / Aqij^ evolution terms are small in DIS, even for ~ 2 GeV^ [2^, 
are well on our way to proving the dominance of the leading twist term of the amplitudes 
( "handbag approximation") where the virtual photon scatters off a single parton. This 
approximation is a corner stone of the study of the nucleon structure in terms of GPDs. 



C. Physics Goals and Proposed DVCS Kinematics 

The present proposal cannot fully disentangle the spin-fiavor structure of the GPDs. We 
itemize here the measurements we will perform and the physical insights we expect to obtain. 

• Measure the ep epry cross sections at fixed xbj over as wide a range in as possible 
for A; < 11 GeV. This will determine with what precision the handbag amplitude 
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dominates (or not) over the higher twist amplitudes. More generally, we consider the 
virtual photon at high as a superposition of a point-like elementary photon and a 
'hadronic' photon (gg, vector mesons) with a typical hadronic transverse size. The 
dependence of the DVCS cross sections measures the relative importance of these two 
components of the photon 

• Extract all kinematically independent observables (unpolarized target) for each Q^, 
XBj, t point. These observables are the angular harmonic superposition of Comp- 
ton Form Factors (CFFs). As functions of 0^^, the observable terms are cos(n0^^) 
for n G {0,1,2}, and sin(?7,0^^) for n G {1,2}, with additional 1/[J + -R'cos(</)^^)] 
modulations from the electron propagators of the BH amplitude ( 311 B|) . Each of 
these five observables isolates the 3fJe or 3m parts of a distinct combination of linear 
{BH ■ DVCS^) and bilinear {DVCS ■ DVCS^) terms. 

• Measure the t-dependence of each angular harmonic term. The t-dependence of each 
CFF is Fourier-conjugate to the spatial distribution of the corresponding superposition 
of quark distributions in the nucleon, as a function of quark light-cone momentum- 
fraction. In a single experiment, we cannot access this Fourier transform directly, 
because we measure a superposition of terms. However, we still expect to observe 
changes in the t-dependence of our observables as a function of xbj ■ In particular, the 
r.m.s. impact parameter of a quark of momentum fraction x must diminish as a power 
of (1 — x) as a; — > 1. This is not a small effect, between x = 0.36 and x = 0.6, we 
expect a change in slope (as a function of t) of a factor two in individual GPDs. 

• Measure the ep — ^ epn^ cross section in the same kinematics as DVCS. 

This will test the factorization dominance of meson electro-production. The longi- 
tudinal cross section {dPai) is the only leading twist (twist-2) term in the electro- 
production cross section. In this experiment, we do not propose Rosenbluth separa- 
tions of (fcTL from (Par- However, as a function of , the ratio (fax/d'^o-L falls at 
least oc l/Q^- Thus the handbag contribution to cPcl can be extracted, within statis- 
tical errors, from a l/Q"^ expansion. The air, o-tt, and (Tlt' terms will be obtained 
from a Fourier decomposition of the azimuthal dependence of the cross section. These 
observables will provide additional constraints on both the longitudinal and transverse 
currents in pion electro-production. 
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DVCS measurements in Hall A/JLab 




Unphysical with E < 11 GeV 

beam 

E,_= 6.6 GeV 



E,,,,= 8.8 GeV 



] E,_= 11.0 GeV 
= 5.75 GeV 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9, 



Figure 4: Proposed DVCS kinematics for H(e, e'^)p measurements in Hall A with 3, 4, and 5 pass 
beams of CEBAF at 12 GeV. The diamond shapes trace the approximate acceptance of the HRS 
in each setting. The boundary of the unphysical region corresponds to the maximum possible 
at a given xbj for 11 GeV. This corresponds to 180° electron scattering, equivalent to 9g = 0° . The 
points at E-Qeam = 5.75 GeV were obtained in EOO-110. 



The handbag amplitude of pion electro-production is the convolution of the axial 
GPDs H and E with the pion distribution amplitude (DA) For charged pion 
electro-production, the E term is expected to be dominated by the pion-pole-this is 
the mechanism used to measure the pion form factor in electro-production on the pro- 
ton. However, neutral pion electro-production will be dominated by the non-pion-pole 
contributions to H and E. 
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To study the DVCS process with CEBAF at 12 GeV, we limit our definition of DVCS to 
the following kinematic range for the central values of our Hall A configurations: 

Q2 > 2GeV2 



Our proposed kinematics, and the physics constraints of Eq. [T]are illustrated in Fig. HI We 
propose a scan of the cross sections for three values of a;^- Our maximum central of 
9 GeV^ is higher than the range of the t-distributions of published DVCS cross sections 
(integrated over 0^^) from HERA for xbj ~ 10^^. The beam time estimates for these high 
statistics measurements are listed in Table HI 

Within the experimental constraints which we detail in section IIIH our choice of kine- 
matics responds to the physics goals of high beam energy DVCS measurements. 

• At each xbj point, we measure the maximum possible range in Q^. Although our 
preliminary data from EOO-110 at Xbj = 0.36 show no indication of dependence in 
the observable 53m[C-^(jF)], this could result from a compensation of terms of higher 
twist and QCD evolution. In this proposal, we will double the range in at xbj = 
0.36, and provide a nearly equal range at xbj = 0.5 and 0.6. 

• Simple kinematics dictates that as the momentum fraction x + ^ of the struck quark 
goes to 1, its impact parameter b/{l + ^) relative to the center of momentum of the 
initial proton must shrink towards zero. Burkardt has discussed the "shrinkage" of the 
proton as X ^ 1 in the limit ^ = 0. j^, The perpendicular momentum transfer 

is Fourier conjugate to the impact parameter b. As a function of x, the mean 
square transverse separation between the struck quark and the center-of-momentum 
of the spectator system is 



where (6^) is the mean square impact parameter of the struck quark. Burkardt con- 
sidered GPD models of the form 



-t < 



(1) 



(rl)(x) = (6^)/(l-x)^ 



(2) 



Hf{x, 0, t) = qf{x) exp [at(l - 



(3) 
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k 


Beam Time 


Total (Days) 


(GeV2) 


(GeV) 


XBi = 0.36 xb\ = 0.50 


XBi = 0.6 (Days) 


3.0 


6.6 


3 




4.0 


8.8 


2 




4.55 


11.0 


1 




3.1 


6.6 


5 




4.8 


8.8 


4 




6.3 


11.0 


4 




7.2 


11.0 


7 




5.1 


8.8 




13 


6.0 


8.8 




16 


7.7 


11.0 




13 


9.0 


11.0 




20 


TOTAL 




6 20 


62 88 



Table I: Proposed Beam Time as a function of {Q^, k, xbj) DVCS kinematics. 

The bound that (r^) remain finite as x 1 requires n> 2. On the other hand, n = 1 
has been used extensively for modeling GPDs. At x = 0.6, the choice of n = 1 or 
n = 2 changes the logarithmic slope d In H/dt by a factor 2.5. This model illustrates 
that without measurements, there are very large uncertainties in the behavior of the 
GPDs at large x (independent of the t = constraints), and that measurements will 
improve our understanding of the transverse distance scales of the quarks and gluons 
inside the proton. 

Our specific choice of kinematics in Fig. HI and Table HI maximizes our range in xbj , while 
maintaining measurements as a function of for each value of xbj . 
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II. 



DVCS OBSERVABLES 



A. DVCS Kinematics and Definitions 

Fig. m defines our kinematic four- vectors for the ep — *■ epj reaction. The cross section is 
a function of the following invariants defined by the electron scattering kinematics: 



y 



{k + pf = M'^ + 2kM 
{q + pf = M^ + Q^ 



1 



- 1 



2q.p W'^-M'^ + Q'^ 
q ■ p/k ■ p. 



(4) 



The DVCS cross section also depends on variables specific to deeply virtual electro- 
production: the DVCS scaling variable ^, the invariant momentum transfer squared t to 
the proton, azimuth of the final photon around the q-vector direction. The DVCS scal- 
ing variable ^ is defined in terms of the symmetrized momenta: 

^2 



q-P 



t/Q' 



■J2-XBj(i-t/g2) 



q^ 



2 - XBj 

(g + g')72 

{p + pT- 

AM^ -t 
{p - pf 



for |t|/Q2 « 1 



(5) 



(6) 



In light cone coordinates {p^-,'P±-,p ), with p^ = ±p^)/v^, the four-vectors are: 



p 

P' 
P 



l^ + U^ , 2 ' 2P+(l+0 
, 2 ' 2P+(l-0 



P+,0, 



4M2 - f 
2P+ 



(7) 
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In Fig. d the initial (final) quark light cone momentum fraction is x ± ^ of the symmetrized 
momentum P and (x ± ± of the initial (final) proton momentum {p'^)- 

Our convention for is defined as the azimuthal angle in a event-by-event spherical- 
polar coordinate system (in the lab): 



= q/|q| 

Vq = [kxkVlkxk'l 

Xq = yqy< Zq. (8) 

tan(0^^) = (q' ■ jjq)/ (q' • Xq) 

= [|q|q'-(kxk']/[qxq')-(kxk')]. (9) 

The quadrant of is defined by the signs of (q' ■ yq) and (q' ■ Xq) We also utilize the 
laboratory angle 9.y^ between the q-vector and q'-directions: 

cos6'^^ = ig-q'/|q'|. (10) 

The impact parameter b of the light-cone matrix element is Fourier conjugate to the A^, 
the momentum transfer perpendicular to the light cone direction defined by P^: 



(1-^2) 

(t — tmin) _^ (■'"-'-) 



tmin = t{e,, = 0°) = ^_^2 - [i-XBiii-Myq^)]- ^^^^ 

B. DVCS Cross Section 

The following equations reproduce the consistent expansion of the DVCS cross section 
to order twist-3 of Belitsky, Miiller, and Kirchner lal. Note that our definition of d^-y agrees 

n n n 

with the "Trento-Convention" for (f) [23|, and is the definition used in and Note also 
that this azimuth convention differs from defined in by (p-y-y = iT — (f^. In the following 
expressions, we utilize the differential phase space element d^^ = dQ'^dxBjd(f)edtd(f)^y. The 
helicity-dependent (A = ±1) cross section for a lepton of charge ±e on an unpolarized target 
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is: 



d^a{X,±e 



dcTQ 



dan 



DVCS/ 



dQ'^dxBj 
dao 



dQ'^dx 



Bj L 



J 



(13) 



a: 



QED 



dQ'^dx 



Bj 



The |T'^^p term is given in 1^, Eq. 25, and will not be reproduced here, except to note 
that it depends on bilinear combinations of the ordinary elastic form factors Fi(t) and F2(t). 
The interference term X is: 

3 



Cq + ^^( — 1)" [c^l-^) cos(r;,0^^) — Xs^sm{n(j)^-y 



n=l 



The (—1)" and (—A) factors are introduced by our convention for 0^^, 
bilinear DVCS terms have a similar form: 

2 



relative to 




It 



DVCS 



(A)r 



1 



c^vcs ^ ^ Asr^^sin( 



(16) 



n=l 



The Fourier coefficients Cn and s„ will be defined below (Eq. I^TH?7jl . The and are 
linear in the GPDs, the c^^'"^ and s^^^^ are bi-linear in the GPDs (and their higher twist 
extensions). All of the 0^^-dependence of the cross section is now explicit in Eq. IT5l and fT6l 

The Vi{(f)^.y) are the electron propagators of the BH amplitude: 

Q'n(077) = (k-q'f 
{k' + q'f. 



Q'^2(077) 



(17) 



After some algebra, the kinematic dependence of the pre-factor of Eq. El is more apparent: 



1 



2\2/ 



M2)/t 



XBiy^Vi{(j)^^)V2{(j)^^)t g2 (J _ 2K cos0^^) (1 + J + (t/g2) _ 2K cos< 



(18) 



J 



21 



ye 



1 — w 

^ 2 



1 + 7^ +(l-^Bj)(2-y) 



1 - a^Bj) 
' t 



t 

as t/Q^ 



-t 



VTT 



l-XBj)[l-y] 



(19) 

XBjW^ (tmin - t) 
- M2 Q2 

(20) 
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The coefficient K appears not only in the BH propagators, but also in the kinematic pref- 
actors of the Fourier decomposition of the cross section (see below). For fixed and xbj, 
K depends on t as (tmin — t)/Q'^- 



C. Fourier Coefficients and Angular Harmonics 



The Fourier coefficients c£ and of the interference terms are: 

t 



4 



8(2 - y)^e { ^1-J^K'C^{J^) + ^(1 - y){l - xbj) [C^ + AC^] (^) 



-8K 



y 

(2-2y + y2) 
-Ay(2-y) 
-IQK^ f (2 - y) W 3?e 
-\y J \ $5m 



3?e 



$5m 



C^(^) 



(21) 



2-a;Bj 

The Fourier coefficients cf , sf are gluon transversity terms. We expect these to be very small 
in our kinematics, though it would be exciting if they generated a measureable signal. The 
and AC^ amplitudes are the angular harmonic terms defined in Eqs. 69 and 72 of (we 
have suppressed the subscript "unp" since our measurements are only with an unpolarized 
target). These angular harmonics depend on the interference of the BH amplitude with the 
set JF = {H, £, H, S} of twist-2 Compton form factors (CFFs) or the related set J-''^^ of 
effective twist-3 CFFs: 



t 



4M2 



(22) 
(23) 



[C' + AC^] (^) = F,{t)n{^, t) - -^F2{t)S{^, t) - eCMit) im, t) + 8{^, t)] . (24) 

Note that [C^ + AC^] depends only on H and S. The usual proton elastic form factors, 
Fi, F2 and Gm = Fi + F2 are defined to have negative arguments in the space-like regime. 
The Compton form factors are defined in terms of the_ vector GPDs H and E, and the axial 
vector GPDs H and E. For example (/ G {u, d, s}) 



E 



e 



"+1 

+V I dx 



m [Hf{^,^,t)-Hj{-^,U)] 



(25) 
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Twist-3 CFFs contain Wandzura-Wilzcek terms, determined by the twist-2 matrix elements, 
and dynamic qGq twist-3 matrix elements. The twist-2 and twist-3 CFFs are matrix elements 
of quark-gluon operators and are independent of (up to logarithmic QCD evolution). The 
kinematic suppression of the twist-3 (and higher) terms is expressed in powers of —t/Q"^ and 
(^min — t)/Q'^ in G..g. the i^-factor. This kinematic suppression is also a consequence of the 
fact that the the twist-3 terms couple to the longitudinal polarization of the virtual photon. 
The bilinear DVCS Fourier coefficients are: 

cr'''' = 2(2-2y + y2)cDVCS(_^^_^*) 

„DVCS 



-.T>VCS 



SK 




II 


' 3fJe 1 


2-XBi 






3m 1 



The cP^*""^ coefficient is a gluon transversity term. 
The DVCS angular harmonics are 

cDvcs^jr^ ^*) = _^ U{i - XBj) (nn* + nn*) - x%2^e 



C^VCS^jreff^jr*^ (26) 



ns* + ns* 



- {^h + (2 - ^B.)^4^) SS* - • (27) 

The twist-3 term C^^^^(J^"^, J^*) has an identical form, with one CFF factor replaced with 
the set J^"^^. The C^^'~^^{J^t,^*), appearing with a cos(20^^) weighting, also has the same 
form as Eq.[23 but now with one set JF replaced by the set JF-p of (twist-2) gluon transversity 
Compton form factors. 



D. BH DVCS Interference and Bilinear DVCS Terms 

The BH-DVCS interference terms are not fully separable from the bilinear DVCS terms. 

) in the general form 

+AJCn53[C^''=^P]sin(n0^^)} (28) 

The factors Kcn,sn are the purely kinematic pre-factors defined in Eqs. [I3H2I1 The experi- 
mental coefficients 3?e, 3mC,^'°^P include contributions from the bilinear DVCS terms, that 
mix into different orders in cos(n0) or sin(n0) due to the absence of the BH propagators 
V1P2 in the DVCS^ cross section. 
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From Eq. CHI [13 and we obtain the generic enhancement of the interference terms 
over the DVCS^DVCS terms (of the same order in sin^ncj)^^) or cos{n(f).y^)): 

BH-DVCS 



oc ^ -■ (29) 



DVCSWVCS 

In each setting (xbj, Q^) setting, for each bin in t, we therefore have the following exper- 
imental Twist-2 DVCS observables 

53m[C^'^^P(^)] = Sm[C^(^)] + (r/,i)!3m[C^^^^(^*,^°fif)] (30) 
3?e {[C + AC]^'^^P(^)} = 3fJe {[C^ + AC^](^)} + (r/o) 3?e {C^^'^'^{T\T)\ (31) 
sfte {C^'^"P(^)} = 3?e {C^(^)} + (r/ei) 3?e {C°v^S(^*, ^)} . (32) 

The coefficients (?7a) are the acceptance averaged ratios of the kinematic coefficients of the 
bilinear DVCS terms to the BH-DVCS terms. In addition, we have the experimental Twist-3 
DVCS observables: 

Sm[C^'''"P(J-^^)] = 3m[C^(^°^)] + (r/,2)»m[C^^^^(^*,^^'^))] (33) 
3?e[C^'^"P(J^^^)] = 3?e[C^(J^^^)] + (r/,2)3?e[C^^^^(J^*,J^^'^)]. (34) 

The values of the ?7a coefficients in the EOO-110 kinematics are summarized in Table ITTl 
They are small, though they grow with \t\. The bilinear term in Eq. [HOlis a Twist-3 observ- 
able, therefore the coefficient (r/^i) will decrease as 1/ a/Q^. Based on the values of Table ITU 
and using our results in EOO-110 to estimate 53m[C^^'"'^(jF*, JF'^^))], we conclude that the 
bilinear term likely makes less than a 10% contribution to 3m[C'^*'^'^(JF)]. In any case, any 
comparison of the experimental results with model calculations, or fit of model GPDs to the 
observables, must include the bilinear terms, with the experimental values of (?7a). 



III. DESCRIPTION OF EXPERIMENT APPARATUS 

This proposal is based directly on the experience of EOO-110. We present a sketch of the 
DVCS layout in Hall A in Fig. El We use the standard 15 cm liquid hydrogen target. We 
detect the electrons in the HRS-L and photons (and 77) in a PbF2 calorimeter at 

beam right. We note in Fig. El the modified scattering chamber from EOO-110 and a new 
modified downstream beam pipe. The scattering chamber is 63 cm in radius, with a 1 cm Al 
spherical wall facing the PbF2 calorimeter and a thin window (16 mil Al) facing the HRS-L. 
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The HRS-Left (HRS-L) limits the central values of the scattered electron momentum to 

k' < 4.3GeV/c, 

Oe > 12.5°. (35) 

As detailed in section lTlI (][ to handle both the instantaneous pile-up and integrated radiation 
dose in the calorimeter, we limit the placement of the calorimeter to 

> r. (36) 

At the same time, to ensure adequate azimuthal coverage in (j)^^, we limit the placement of 
the calorimeter 

Og > 10°. (37) 

The following subsections detail our technical solutions, and demonstrate that these tech- 
nical constraints do not limit the physics scope of this proposal. 

A. High Resolution Spectrometer 

Our proposed kinematics, and the physics constraints of Eq. [T] are illustrated in Fig. |3| 
The individual beam energies are illustrated in Fig. El with the HRS constraints (Eq. ESI) 
superimposed. We note the following points with regard to the HRS and calorimeter (Eq.|H3) 
constraints: 

• The HRS constraints have no effect for > 2 GeV^ at k = 6.6 GeV. 



Table II: Weighting factors of bilinear DVCS terms for BH-DVCS observables in EOO-110. 
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at 3.00 M 

Figure 5: Hall A layout for DVCS with CEBAF at 12 GeV. The scattering chamber is identical to 
the EOO-110 chamber, with a midplane 63 cm radius spherical section with 1 cm Al wall thickness, 
and a 16 mil Al window facing the HRS L. We propose a conical downstream beam pipe, with 
half-opening angle of 6 degrees on beam-right, and length 3m. The drawing shows the expanded 
PbF2 calorimeter at beam-right, in the closest and farthest configurations: front face 150 cm and 
300 cm, respectively, from target center. The calorimeter is shown in its smallest angle setting 
(inner edge at 7°). 

• At k = 8.8 and 11 GeV, primarily the k < 4.3 GeV constraint removes roughly the 
lovi^er half of the range at each xbj . This range for each xbj is covered by the 
lower beam energies. 

• According to the Base Equipment plan for CEBAF, the three Halls A, B, C, when 
running simultaneously, must operate at different multiples of 2.2 GeV. 

• The small angle calorimeter cut of 6*^ > 10° removes the very highest physically allowed 

values at each xbj. 



21 



Table III: Detailed DVCS Kinematics. The first line is from EOO-110, and is included for comparison 
purposes. The angle 9q is the central angle of the virtual photon direction q= {k — k'). 
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• All of the constraints prevent us from kinematics for xbj < 0.2. 

• Kinematics at xbj = 0.7 are allowed at k = 8.8 and 11 GeV. At this time, it is very 
difficult to make reliable estimates of the DVCS signal at this extreme xbj, and we do 
not include these kinematics in our proposal. 

The detailed kinematics are summarized in Table ITTTl 

The HRS performance is central to this experiment. Once we make a selection of ex- 
clusive H(e, e''j)p events, the resolution in t = (g — q'Y is determined by approximately 
equal contributions from the HRS momentum resolution and the angular precision of the 
direction q'. The resolution on the direction of q' has equal contributions from the position 
resolution in the calorimeter and the vertex resolution, as obtained from the spectrometer. 
The combination of precise acceptance, kinematic, and vertex resolution in the HRS makes 
it possible to make precision extractions of the DVCS observables. 
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Figure 6: Proposed DVCS central kinematics for H(e, e'j)p measurements in Hall A with 6.6, 
8.8, and 11 GeV incident beam. The experimental constraints of Eg. I35|371 are indicated. Each 
experimental point is the center of one of the 'diamond' regions of Fig. ^ 



B. Beam Line 



In EOO-110 and E03-106, with the front face of the calorimeter 1.10 m from target center, 
we had a 6" diameter cylindrical beam pipe, welded directly to the scattering chamber. This 
aperture corresponds to roughly a 6° half-opening angle from target center. We require 
placing the calorimeter at distances from 1.5 m to 3.0 m from the target f ^III Cj) . To avoid 
excess calorimeter background generated by secondaries striking the beam pipe, we propose 
a thin walled (1/8" AI) conical beam pipe, welded to the same aperture. The cone must 
be slightly off-axis, to prevent an interference with HRS-Ql at the minimum spectrometer 
angle setting of 18° (Table ITTT| . while at the same time preserving a full 6° aperture on 
Beam-right. The cone should be 3.0 m long, to continue downstream of the calorimeter. 
This will require moving or replacing a gate valve and pumping station. 
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C. PbF2 Calorimeter 



Table IV: Properties of PbF2. 



Density 7.77 g/cm'^ 

Radiation Length 0.93 cm 

Moliere Radius 2.20 cm 

Index of Refraction (A = 180 nm) 2.05 

(A = 400 nm) 1.82 

Critical Energy 9.04 MeV 



We will detect the scattered photon in a 13 x 16 element PbF2 calorimeter. This is the 
existing 11 x 12 EOO-110 calorimeter, with 76 additional elements. Each block is 3 x 3 
cm^ X 20Xq. The additional blocks will add two more rows on the top and bottom, and two 
columns on the wide angle side. The properties of PbF2 are summarized in Table Hvl The 
important design considerations for DVCS are as follows. 



• PbF2 is a radiation hard pure Cerenkov crystal medium 




• With no scintillation light the calorimeter signal is insensitive to low energy nuclear 
particles, and the pulse rise and fall time is determined only by geometry and the 
response of the PMT. This allows us to use the IGHz Analog Ring Sampler (ARS) 
digitizerj^^ to minimize pileup. 

• The high luminosity of this proposal requires fast response PMTs operated at low 
gain and capacitively coupled to a pre-amplifier. The low gain reduces the DC anode 
current. The capacitive coupling removes the average pile-up from low energy 7-rays. 

• The small Moliere radius (2.2 cm) allows us to separate closely space showers from 7r° 
decay, and minimize shower leakage at the boundary. 

• The short radiation length minimizes fluctuations in light collection from fluctuations 
in the longitudinal profile of the shower. 

• The low value 9 MeV of the critical energy (roughly the energy threshold for which 
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bremsstrahlung energy loss exceeds ionization loss for electrons) also improves the 
resolution e.g. relative to Pb-Glass. 

• In EOO-110, we obtained a signal of 1 photo-electron per MeV of deposited energy in 
the E.M. shower, and a energy resolution of 2.4% from elastic H(e, e'^aioPuRs) elec- 
tron of 4.2 GeV. For our simulations, we project a resolution oi ce/E = 2.0% © 
(3.2%)a/ (1 GeV)/q'. We also achieved a spatial resolution of 2 mm at 4.2 GeV. From 
the combination of energy and spatial resolution, we obtained a 7r° — 77 mass reso- 
lution of 9 MeV. 

The size, granularity, and position of the calorimeter must accommodate the following con- 
straints: 

1. Nearly 27i azimuthal photon acceptance for \A±\ < 0.6 GeV/c. independent of the 
central kinematics. The angular size required therefore shrinks as A™'^^/g'(0°). 

2. Good separation of the two clusters from the 7r° — 77 decay, in order to measure 
the H(e, eV°)p reaction. For tt" 77 reconstruction, we require a center to center 
separation of 3 PbF2 blocks, or 9 cm. In high energy DVCS or deep virtual 7r° kine- 
matics, q'dycs ~ -^T- In this limit, the minimum half opening angle of the vr^ — s> 77 
decay is 6'^^ > niT^/q'. The minimum distance of the calorimeter from the target, 
based on the cluster separation requirement is given as the D parameter in Table d 
Furthermore, 50% of the 7r° decay events yield both photons within a half angle cone 
of < {^/Zj2){m^/q'). 

3. Maximum distance from target to minimize pile-up and radiation dose per block at 
fixed luminosity. 

Item [2I requires us to increase D in proportion to q' . However, from item Q] we see that 
the acceptance is remains invariant. The solid angle per PbF2 block at the distance D 
determines the maximum feasible luminosity for each setting. Therefore item IHl allows us to 
increase the luminosity in proportion to oc g'^. 

We will calibrate the central blocks of the calorimeter via elastic H(e, (^'caioP'RRs) mea- 
surements We anticipate three sequences of elastic measurements of one day each at the 
beginning, middle, and end of each scheduled period of running time. We will then cross 
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calibrate all of the blocks and maintain a continuous monitor of the calibration with the vr 
mass reconstruction from H(e, e'n^)X events. We anticipate sufficient statistics to obtain 
an independent calibration from each day of running. The elastic calibrations also serve to 
verify the geometrical surveys of the spectrometer and calorimeter. 

D. Luminosity Limits and Optical Curing of Calorimeter 

In E03-106, we took D{e,e''y)X data at a maximum luminosity (per nucleon) of 4 • 
10^^/cm^/s, with the front face of the PbF2 calorimeter at 110 cm from the target center. 
We did not obtain any degradation of resolution in the calorimeter from pileup of signals. 
During the entire 80 day run of EOO-110 and E03-106, we delivered a total of 12 C to the 
15 cm liquid Hydrogen and Deuterium targets. We performed absolute calibrations of the 
calorimeter with elastic H(e, e'^^i^pim^s) events at the beginning and end of data taking. We 
observed up to 20% decrease in signal amplitude in individual blocks, without observable 
loss in missing mass resolution after recalibration. Custom pre-amplifiers were used to keep 
the PMT anode current small. Therefore the loss of amplitude is attributed to degradation 
of the transmission properties of the blocks, and not to degradation of the photo-cathodes 
of the PMTs. Independent numerical and analytic simulations indicate that below 10 
degrees, the radiation dose to the calorimeter is dominated by M0ller electrons (and related 
bremsstrahlung). Beyond 20 degrees, the dominant background arises from decay photons 
from inclusive 7r° photo-production. The simulations also indicate that from 7.5 to 11.5 
degrees, the radiation dose diminishes by a factor 5. On this basis, we conclude that roughly 
50% of the radiation dose received by the small angle blocks of the calorimeter occurred when 
the small angle edge of the calorimeter was at 7°. This is the same minimum angle we will 
use in the present proposal. 

Achenbach et al. studied the radiation damage and optical curing of PbF2 crystals 
They found the radiation damage to be linear for doses up to 8 KGy (from a ^°Co source). 
For a dose of 1 KGy, they observed a loss of 25% in transmission for blue light of A = 400 
nm. They also obtained good results for curing the radiation damage by exposure to blue 
light. The front face of the 1000 element A4 array was exposed for 17 hours to a Hg(Ar) 
pencil lamp (filtered to pass only A > 365 nm) at a distance of 50 cm (intensity on the 
calorimeter surface of 2/iW/cm^). The transmission of 400 nm and 330 nm light returned 
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to 100% and 97%, respectively, of its initial value. 

At this stage, it is not possible to compare the absolute dose in our simulations with 
the dose recorded in the MAMI-A4 trials. The A4 dose is recorded as a volumetric dose 
(IGy = IJoule/kg) yet the gamma rays from ^°Co are predominantly absorbed in a layer 
of thickness 1/10 the transverse size of the crystals. The radiation dose in Hall A during a 
DVCS experiment is primarily from photons and electrons 50-1000 MeV. Thus we consider 
the absolute scale of dose comparison to be uncertain by a factor of 10. Instead, we normalize 
the future radiation damage of the calorimeter to the maximum value for the fractional signal 
attenuation per integrated luminosity obtained in conditions during EOO-110 and E03-106 
that were nearly identical our proposed configuration. We utilize the studies from MAMI-A4 
to determine that a radiation-dose induced attenuation of up to 25% (A = 400 nm) can be 
cured (to within 1%) with a 17 hour exposure to blue light. This curing must be followed 
by several dark hours to allow the phosphorescence of the PMT photo-cathodes to decay. 

Pile-up within the 20 ns analysis window of the pulse shape analysis of the PbF2 signals 
will limit our instantaneous luminosity. Based on our previous experience we can operate 
at an instantaneous luminosity times acceptance per PbF2 block of 
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110 cm 



(38) 



with D the distance from target to calorimeter. Our projected count rates and beam times 
are based on this luminosity as a function of the kinematic setting. 

We plan to use blue light curing of the blocks every time the signal attenuation reaches 
20%. For those settings with the minimum edge of the calorimeter 6q]^^ equal to 7°, this 
corresponds to 5 days of running at the luminosity of Eq. EHl We will use the 7r° — > 77 mass 
resolution from both single arm H(e, 7r°)F (prescaled) and coincidence H(e, e'7r°)X events 
to monitor the light yields in the PbF2 array. At larger calorimeter angle settings, the time 
between curing will be correspondingly longer. We estimate a total of 10-12 curing days 
during the experiment (in addition to the running time). 

E. Trigger 



The electron detector stack in the HRS will be the standard configuration of VDC I and 
II, segmented SI and S2, gas Cerenkov, and Pb-Glass calorimeter. The Cerenkov {Cer) and 
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Pb-Glass provide redundant electron identification in the ofi'-line analysis. The main HRS 
trigger is 

HRS = [U,j{SlRinSlLi)n{S2Rjr]S2Lj)]r]Cer. (39) 

This requires a coincidence between the PMTs at the two ends of at least one scintillator 
paddle in each of 5*1 and 5*2. In addition, we require a coincidence with the Cerenkov 
counter. Supplementary prescaled triggers with either 51, 5*2, or Cer removed from the 
coincidence will monitor the efficiency of each trigger detector. 

In addition to the HRS electron trigger, we will upgrade our present coincidence 
validation/fast-clear logic. The calorimeter signals are continuously recorded by the 128 
sample® 1 GHz ARS 0Q 

array. This array is digitized, at a slower rate, following a trig- 
ger validation. For each HRS trigger, we stop the ARS sampling, and trigger a Sample and 
Hold (SH) circuit, that is coupled to the calorimeter signals via a high impedance input (in 
the future, we may replace this with a pipeline ADC). The SH signals are digitized and then 
in a field programmable gate array (FPGA), we will form the following trigger validation 
signals. 

1. A validation of the HRS-electron trigger, based on detecting at least one 2x2 cluster 
above a programmable threshold E"^^. 

2. A 7r° trigger based on detecting at least two separated 2x2 clusters, with each cluster 
above a programmable threshold Ei and the sum of the two clusters above a pro- 
grammable threshold E^,. This trigger will select candidate H(e, eV°)X events. 

If a valid signal is found, the ARS array is digitized and recorded (together with the HRS 
signals). If the HRS trigger is not validated by the calorimeter signals, a fast clear is 
issued to the ARS array, no digitization occurs, and acquisition resumes. During EOO-110, 
the SH/Fast-clear cycle took 500 ns. With upgrades to the FPGA, we can shorten this 
deadtime by a factor of two. In addition, upgrades to the VME standard will allow us to 
increase the total bandwidth of data acquisition for this proposal. 



IV. PROJECTED RESULTS 



The detailed (e, e') kinematics, calorimeter configuration, H(e, e'j)X missing mass reso- 
lution in M|^, count-rates, and beam time are summarized in Table Ivl The beam time is 
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Table V: Experimental Conditions for DVCS. For each (Q^,fc,rEBj) setting, we present: Maximum 
photon energy q'{0°); Calorimeter distance D] Virtual photon direction Oq] Angle 0"^]^ of the edge 
of the calorimeter, relative to the beam line; Kinematic niinimuni |^min| ^^nd upper bound |tmax| 
for approximately full acceptance in (j)^^; Resolution cr(M^) in H(e, e'7)X missing mass squared; 
Luminosity £; H(e, e)X inclusive trigger rate, H(e, e'^)p exclusive DVCS count rate. The beam 
time is calculated to obtain an estimated 250K events at each setting, or at least 40,000 events 
per bin in tmm — t- The distance D of the front face of the calorimeter from the target center 
is optimized for the separation of the clusters from vr^ 77 decay (section IlII Cp . The intrinsic 
missing mass resolution in EOO-110 is cr(Mj^) = 0.20 GeV^. The luminosity is determined by the 



maximum rate allowed by pileup in the calorimeter. This luminosity is proportional to ( 8III Pp . 
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2.5 


13.1 


8.6 


-0.69 


-1.10 


0.20 


2. 


274 


0.28 


13 


9.0 


11.0 


0.6 


7.62 


3.0 


10.2 


7.3 


-0.71 


-1.14 


0.22 


3. 


117 


0.17 


20 



chosen to give a total statistics of 250K events per (xBj,Q^) setting. 



A. Systematic (Instrumental) Errors 

Table ED shows the main systematic errors on the cross-sections extraction during esti- 
mation of those errors for the proposed experiment. The main improvements of the errors 
are due to : 
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- The beam polarization measurements: the upgrade of the Compton polarimeter along 
the Hall A beam line will reduce the relative precision AP/P from 2% to 1%. The 
average beam polarization for EOO-110 was 75% and is expected to be 85% for the 
proposed experiment. 

- Drift chamber multi-tracks : the number of events for which multiple tracks are re- 
constructed in the focal plane of the HRS is proportional to the single rates. Multiple 
tracks events are discarded during the analysis and the live time of the experiment 
corrected accordingly. Because the proposed experiment is measured deeper in the 
DIS region, we expect that the correction will be smaller as well as its uncertainty. 

- TT^ subtraction : the limitation on the precision of the correction for EOO-110 was both 
statistical and limited by the kinematical range of the recorded 7r° events. For the 
proposed experiment, we plan to improve the trigger such that more 7r° events are 
recorded. A special trigger with lower threshold but the requirement of two clusters 
being above it, will be implemented. 

- e(p,e'7)7rN contamination : the missing mass squared resolution is the key ingredient 
for this contamination. The proposed estimation takes into account the variation of 
this resolution for the new kinematics. 

The systematic errors in Table IVII are summarized separately for the cross section sum 
d'a + (io^, and for the cross section difference d'a — d^ . 

B. Statistics 

Fig. EHHl show acceptance, counting statistics, and cross section distributions for specific 
(xBj,Q^) kinematics. In each figure, we present the cross section weighted acceptance dis- 
tribution of a basic set of kinematic variables in the upper left panel. The central plot of 
this panel is the H(e, e'7)X missing mass squared Mj. distribution for simulated exclusive 
events. The blue curve is a simulation for EOO-110. The missing mass resolution a{M]^) for 
each setting is tabulated in TableEl For EOO-110, the resolution is (t(M|) = 0.20 GeV^. In 
the lower right kinematic plot. The green rectangle is the fiducial surface of 11 x 14 blocks 
of the calorimeter (we require the centroid of EM shower to be a a minimum of of one 
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Type 


Relative 


errors (%) 






EOO-110 


proposed 


Luminosity 


target length and beam charge 

o o o 


1 


1 


HRS-Calorimeter 


Drift chamber multi-tracks 


1.5 


1 




Acceptance 


2 


2 




Trigger dead-time 


0.1 


0.1 


DVCS selection 


vr'^ subtraction 


3 


1 




e(p,e'7)7rN contamination 


2 


3 




radiative corrections 


2 


1 


Total cross section sum 


4.9 


4.1 


Beam 


Polarization AP/P 


2 


1 


Total cross section difference 


5.3 


4.2 



Table VI: Relative systematic error budget for EOO-110 and for the proposed experiment. 



block from the calorimeter edge). The points in red on this plot show the distribution of 
directions of the q- vector. The upper right panel of the figures shows the counting statistics 
(helicity sum and helicity difference) as functions of (p-y-y for five bins in t. As —t increases 
from left to right, the simulation shows the gradual loss of acceptance for near 0°, due 
to the off-centered calorimeter. The bottom right panel shows the projected helicity sum 
and helicity difference cross sections, with statistical errors, in the same bins in t and 0^^. 
The statistics are for the projected beam time of this proposal. These estimates were made 
with our simulation code TCHIB. The cross section model for these estimates is described 

in 

With 250K events in each (xBj,Q^) setting, we obtain a high precision determination of 
all the observables detailed in the Section III Dl 
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K=1 1 GeV, Q =9 GeV^x,=0,6, 0,=3O.23°,k'=3 GeV.0„.= -1 1° 
Calo 1 3x16 Blocks at 3 meters ^=2.97X10*" cm"'s'', 400 Hours 
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Figure 7: DVCS Distributions for setting 
A; = 11 GeV, = 9 QeV^, xbj = 0.6. 
Top Left: Cross section weighted acceptance 
distributions. 

Top Right: Helicity sum and hehcity differ- 
ence projected counts as a function of 4)^^ in 
five bins in t. 

Bottom Right: Hehcity sum and hehcity dif- 
ference projected cross sections, with statis- 
tical uncertainties, as functions of (\>^^ in the 
same bins in t. 
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K = 8.S GeV, Q'=4,8 GeV',Xe=0.5, 0,=22.2°,k'=3,S8 GeV.0„,.=- 1 4.47° 
Colo 13x16 Blocks at 2 meters Lu= 1 .32 X 1 0*'' cm"'s'', 85 Hours 
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Figure 8: DVCS Distributions for setting 
k = 8.8 GeV, = 4.8 GeV^, xbj = 0.5. 
Top Left: Cross section weighted acceptance 
distributions. 

Top Right: Helicity sum and hehcity differ- 
ence projected counts as a function of 4>^.y in 
five bins in t. 

Bottom Right: Hehcity sum and hehcity dif- 
ference projected cross sections, with statis- 
tical uncertainties, as functions of (p^ry in the 
same bins in t. 
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K=5,6 GeV, Q'=3,0 GeV',x,=0,36, 0,=25.5°,k'=2.1 5 GeV,0„>,= - 1 4.5° 
Colo 13x16 Blocks at 1.5 meters Lu=0.74X 10*" cm'^s'', 60 Hours 



0.2 0.5 0.8 



0.5 0.75 



0" (GeV") 



10 
(GeV") 



2 
M,' (GeV") 



-2 -1 

t (GeV') 



Gj,"* (degree) 



2.5 5.75 
q' (GeV) 



colorimeter 



-30 30 

Yyx„ (cm/cm) 



Figure 9: DVCS Distributions for setting 
k = 6.6 GeV, = 3.0 GeV^, xbj = 0.36. 
Top Left: Cross section weighted acceptance 
distributions. 

Top Right: Helicity sum and hehcity differ- 
ence projected counts as a function of 4>^.y in 
five bins in t. 

Bottom Right: Hehcity sum and hehcity dif- 
ference projected cross sections, with statis- 
tical uncertainties, as functions of (p^ry in the 
same bins in t. 
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C. Illustration of Prospective Physics 

In each setting (xBj,Q^) setting, we will extract in the t-bins defined in Fig. EHHl the 
experimental Twist-2 and Twist-3 BH-DVCS interference observables defined in Section. FlI Dl 
In addition, we will obtain upper bounds on the gluon transversity terms cos(30^^) and 
sin (30^^). 

We illustrate the possible xej-dependence of the Twist-2 coefficients in Fig. [101 The 
principal-value integrals in the 3?eC-^ and 3?eAC-^ must change sign somewhere in the range 
< ^ < 1. In Fig. El we present the impact on 53m[C'^] of the DGLAP evolution of 
t; g2) in the VGG model 0, Q Q • 

The VGG model uses the double distributions 

for H: 

Ff{(],a,t) = h{f3,a)qf{(3)/f3"'' (41) 
h{P,a) = db) [{1 - - a'Y /{l-lPlf^' (42) 

The limits of integration in Eq. HOl and the normalization constant C(b) are defined in 
. In these illustrative plots, we have taken the profile parameter b = 1 and the Regge 
parameter a' = 0.8 GeV"^. A similar form is used for H, with g/(/3) replaced with Agj(/3). 
We also set = in these figures. For our count rate estimates in the previous section, 
we have used the factorized ansatz which is obtained in the b ^ oo limit. In Fig. [TH 
the DGLAP evolution is applied only to the parton distributions g/(/5). We note that 
the experimental contribution of the Twist-3 term (?7i)55m[CP^*"^] fEq. l30|) will decrease in 
proportion to 1/ \/Q^. At the same time, higher twist contributions in 53m[C-^] will decrease 
in proportion to l/Q"^ (or higher powers). Thus it should be possible to extract the Twist-2 
and Twist-3 contributions from the Q^-dependence, with only a small model dependence 
from the distinct QCD evolution of the three contributing GPDs. 

Fig. US illustrates the t-dependence of 53mC-^(jF), plotted here as a function of tmin — t 
for our three values of xbj. The t-dependence of this observable is the bilinear combination 
of the Compton form factors with the proton form factors in the BH amplitude. 
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Figure 10: The three Twist-2 observables, as functions of xbj at imm — t = 0.3 GeV^. The model 
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Figure 11: QCD evolution of C-^ , as obtained from DGLAP evolution of H{x,(^,t; Q^) in the VGG 
model Q,y,Q 
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Figure 12: Dependence of 9mC-^(^) as a function of tmm — t for the three values of xbj of this 
proposal at constan t Q ^. The curves are obtained from the contributions of H{x, ^, t) and H{x, ^, t) 
in the VGG model 
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V. SUMMARY 

With a high precision spectrometer, and a compact, hermetic high-performance calorime- 
ter in Hall A, we can obtain high statistics absolute measurements of ep ep-f cross sections 
in a wide DVCS kinematic domain with CEBAF at 12 GeV. We expect to obtain precision 
measurements of three Twist-2 DVCS observables and two Twist-3 DVCS observables. The 

dependence will provide stringent tests of the factorization theorem, and quantify the 
contribution of higher twist terms (which can be modelled as the hadronic content of the 
photon). The t-dependence as a function of xbj (or ^) of the observables will provide our 
first study of transverse profile of the proton as a function of quark light-cone momentum 
fraction. We note that in this kinematic region, models such as the VGG model presented 
in §IV CI must be considered uncertain by at least a factor of two. Our separations of the 
3?e and 53m parts of the DVCS observables will provide a calibration for present and future 
measurements of relative asymmetries. Through the principle value integrals, the 3?e part 
provides access to regions of (x, ^)-space that are not directly accessible via the Q'm part of 
the GPDs. The 3m part observables restrict the GPDs to the sum of the points x = 

We will obtain precision measurements of the cl+o't/^l, o'lt, cttt and ctlt' cross sections 
of neutral pion electro-production H(e, e'7r°)p. The factor is the degree of longitudinal 
polarization of the virtual photon.. The transverse cross section ctt is expected to fall faster 
with (by one power of l/Q"^) than the leading twist term in the longitudinal cross section 
daL- Therefore, from the dependence of the H(e,e'7r°)p cross section, we can extract the 
leading twist (handbag) contribution to this process. 

We require 88 days of production running, with an anticipated additional 12 days inter- 
laced for optical curing of the calorimeter. This will provide a major survey of proton DVCS 
in almost the entire kinematic range accessible with CEBAF at 12 GeV. 
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Appendix A: CONTRIBUTIONS TO HALL A EQUIPMENT FOR 11 GEV 

Our committments to the Hall A Base Equipment are detailed in Table ETH The collab- 
oration will also seek additional funds for the upgrade to the calorimeter and beam line. 



Table VII: Contributions of DVCS Collaboration to Hall A Base Equipment 



Institution 


Funding Source 


Project Description 


Equipment 


LPC-Clermont-Ferrand 


IN2P3-CNRS 


Compton Polarimeter: 


$125K 






HRS Trigger Electronics 


$39K 


Technical Manpower for Base Equipment 




(months) 


LPC-Clermont-Ferrand 


IN2P3-CNRS 


Compton Polarimeter: 








Engineer . 


24 






T echnician 


24 






P hysicist . 


12 






HRS Trigger Electronics 








Engineer . 


6 






T echnician 


6 


Old Dominion U. 


DOE 


HRS DAQ Faculty 


4 




(ongoing grant) 


PostDoc 


12 



Appendix B: PREPARATION OF EXTENSIONS 

We are actively preparing a number of extensions of this proposal, which will lead to a 
comprehensive program of measurements exploiting the unique capabilities of Hall A. 

1. Deuterium 

We currently have preliminary quasi-free and coherent DVCS data from Hall A E03- 
106 in the channels D(e, e'7)pn © D(e, e'7)D. These two channels are (partially) separated 
experimentally by the differential recoil of the coherent deuteron. In a plot of M|^, calculated 
for H(e, e'7)p kinematics, the coherent Deuteron peak appears at Mx = +t/2 (Fig. fT3|l . 
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The smearing of the quasi-elastic D{e,e''y)pn events by the internal momentum distribution 
of the np wavefunction of the deuteron is less than our experimental resolution. Thus our 
isolation of the D(e, e''j)pn channel from the inelastic D(e, e''-f)NNTT channel is only slightly 
degraded compared to the identification of the H(e, e''-f)p exclusive channel. 

In a Quasi-Free (QF) model, the neutron DVCS cross section can be obtained from the 
following subtraction: 

(fa{en e''yn) =(fa{eD e'^np) — d^a{ep e'7p) (Bl) 

All observables on the neutron have the same form as for the proton. Using isospin symmetry, 
the neutron observables interchange the role of up and down quarks. Whereas the proton 
cross sections are four times more sensitive to up than down quarks, the reverse is true for 
the neutron. As a consequence, based on the forward parton limits of the GPDs and the 
Form Factor constraints of the first moments of the GPDs, we expect the neutron observable 
to be dominated by the GPD E. 

In Fig. O we show the helicity correlated statistics in three zones in missing mass of the 
'neutron' spectrum of Eq. IbTI We note that the signal obtained in the 'A^A^-coherent' region 
Mjr < 0.6 GeV^ is of opposite sign from the signal in the Quasi-Free region near M]^ = M^. 

We expect to prepare a separate proposal for £)(e,e'7) measurements for the next 12 
GeV PAC. We expect the Deuterium proposal to include a subset of the present kinematics. 
We would propose to run the two experiments concurrently, with alternating sequences of 
Hydrogen and Deuterium data taking, to minimize systematic errors. 

2. Recoil Polarimetry 

A full DVCS program requires proton polarization measurements, and detection of nu- 
clear recoils in coherent A{e, e'jA) DVCS reactions. The single- and double-spin observables 
of proton recoil polarization in ep — > e'p'j and functionally equivalent to the observables for 
polarized targets: ep e'prf (D. Miiller, private communication 2006). We are working on 
a conceptual design for a large acceptance recoil polarimeter. With an analyzer of 7 cm 
C, we can achieve a Figure of Merit > 0.005 (scattering probability times analyzing power 
squared) for Pp > 550 MeV/c. This could measure both longitudinal and transverse recoil 
proton polarization in DVCS. This would operate in an axial magnetic field at high lumi- 
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Figure 13: Top: Missing mass distribution of E03-106 D(e,e'7)X events, after subtracting a 
normalized ensemble of H(e,e'7)X events. The kinematics are = 1.9 GeV^, a^Bj = 0.36, and 
{t) = —0.3 GeV^. The missing mass is calculated for N{e,e''y)N kinematics. Thus the quasi-free 
n{e,e'j)n peak appears at M'^ = and and the coherent D(e,e'7)D peak is broadened by the 
t-acceptance of the events, and appears at M'^ = + 1/2. In the region > + 1/2, there 
can also be a contribution from D(e,e''y)np on high momentum correlations in the initial state 
of the Deuteron (WA^-Coherenf). The curves are a two-gaussian fit, with the relative positions 
constrained to t/2. The bottom three plots show the helicity weighted statistics N{+) — N{—) in 
the three missing mass regions: 0.0 < Mj. < 0.6 GeV^ (left); 0.0 < M|. < 1.15 GeV^ (center); 
and 0.9 < M'^ < 1.15 GeV^ (right). The WA^-coherent' and 'Quasi-Free Neutron' regions have 
opposite asymmetry. 
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nosity in Hall A, in conjunction with the equipment of this proposal. We intend to have a 
proposal ready for PAC 32. (Summer 2007). 



